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ABSTRACT 


Orifice  and  louver  sound  power  spectra  are  investigated,  using  an 
intensity  probe,  at  various  wind  speeds  in  a  low  noise,  semi-anechoic , 
subsonic  wind  tunnel  for  free  stream  velocities  below  50  meters  per 
second.   The  radiated  noise  is  created  by  turbulent  flow  over  various 
orifice  and  louver  geometries  which  are  flushed  mounted  into  the  wall 
of  a  long  duct.   Five* orifice  samples  of  rectangular  shape  and  various 
transverse  dimensions  as  well  as  four  louver  samples  with  multiple 
rectangular  and  circular  orifices  are  tested.   Also  investigated  is 
the  effect  of  the  leading  and  trailing  edge  angle  on  the  radiated 
sound  power.   The  scaling  laws  of  the  excitation  frequencies  and  the 
speed/power  laws  are  presented  for  ratios  of  the  boundary  layer  thick- 
ness to  the  transverse  orifice  dimension  from  1.01  to  4.29.  A 
detailed  theoretical  model  is  developed  for  rectangular  shaped  aper- 
ture orifices  and  louvers  based  on  the  work  by  Ffowcs  Williams  [9]  , 
Nelson  [11],  and  Corcos  [12].    The  model  describes  the  predominant 
variables  which  effect  the  radiated  flow  noise.   To  ascertain  the  val- 
idity of  the  model,  it  is  used  to  collapse  the  experimental  power 
spectra  for  two  orifice  and  one  louver  geometry.  The  model  showed 
excellent  agreement  with  the  actual  sound  power  measurements.   Direc- 
tivity plots  are  also  presented  to  further  describe  the  orifice/louver 
source . 
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1.   INTRODUCTION 

1 . 1   Background 

Turbulent  flow  over  an  orifice  or  louvered  openings  can  generate 
high  levels  of  broadband  noise.   This  phenomenon  has  been  known  for 
some  time  [1]  and  early  investigators  like  Dunham  [2],  Rossiter  [3] 
and  East  [A]  attempted  to  identify  and  control  the  high  vibrational 
levels  and  possibly  the  radiated  sound  levels  caused  by  the  interac- 
tion of  the  flow  over  these  openings.   Recent  work  by  Elder,  Farabee 
and  DeMetz  [5]  concentrated  on  the  area  of  cavity  resonance  where  the 
opening  provides  the  feedback  path  between  the  closed  acoustic  cavity 
and  the  source,  the  turbulent  flow.   When  the  typical  dimension  of  the 
cavity  is  large  compared  to  the  acoustic  wavelength,  the  source  can 
excite  various  room  modes  of  the  cavity.   If  the  turbulent  vortex 
shedding  frequency  coincides  with  one  of  these  room  modes,  then  a 
strong  feedback  will  occur,  resulting  in  an  acoustic  resonance  of  that 
mode. 

If  the  cavity's  volume  is  further  increased,  multiple  room  modes 
will  be  excited  and  the  strong  feedback  mechanism  will  be  lost.   This 
results  in  a  much  broader  acoustic  response  and  at  lower  amplitude 
levels  than  the  resonant  case  which  is  typically  characterized  by 
large  amplitudes  at  very  discrete  tonal  frequencies.   In  practice,  the 
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detection  of  resonating  cavities  Is  often  easy  due  to  their  strong 
discrete  tones.   However,  once  these  resonances  have  been  corrected, 
or  presumably  prevented  in  the  design,  identifying  the  cause  of  the 
remaining  flow  noise  becomes  much  more  difficult.   Possible  contribu- 
tors to  the  broadband  flow  noise  are:  (1)  the  turbulent  boundary  layer 
direct  transfer  of  energy  to  the  structure  of  the  enclosed  acoustic 
space,  (2)  the  power  flow  through  the  opening  (open  window  effect) 
caused  by  the  acoustic  field  created  by  the  upstream  and  downstream 
turbulent  sources  In  the  boundary  layer,  (3)  the  power  flow  into  the 
acoustic  space  caused  by  the  source  created  in  the  opening  due  to  the 
driving  aerodynamic  field  and  the  shear- layer  edge  interaction  with 
the  downstream  edge  of  the  orifice  and  (4)  the  vibrational  excitation 
of  the  internal  structure  caused  by  the  orifice  source.   The  amount  of 
noise  which  the  latter  contributes  depends  on  the  details  of  the 
internal  structure  and  the  orifice/louver  attachment  design.   There- 
fore to  predict  the  contribution  from  internal  structure  vibration, 
one  must  first  determine  the  orifice  source. 

1.2  Content 

The  contribution  of  the  first  three  possible  contributors  to  the 
total  power  radiated  into  a  closed  acoustic  space  is  analyzed  based  on 
the  results  from  several  experiments.   It  is  shown  that  indeed  the 
orifice  source  provides  both  the  spectral  character  and  the  predom- 
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inant  amount  of  the  total  sound  power.   By  analyzing  and  comparing  the> 
results  with  a  derived  equivalent  source  model,  the  nature  of  the 
source  and  the  parameters  which  affect  its  strength  and  spectral  char- 
acter are  identified.  The  analysis  comes  directly  from  experimental 
measurements  taken  in  a  subsonic  wind  tunnel  at  the  M.I.T.  Acoustic 
and  Vibration  Laboratory. 
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2.   SOUND  INTENSITY  MEASUREMENTS 

2.1  Why  Measure  Sound  Intensity? 

To  describe  an  acoustic  source,  one  must  measure  a  quantity  which 
is  not  dependent  on  the  environment  or  the  surroundings.   The  quantity 
most  often  identified  to  describe  a  source  is  its  power.   Knowing  the 
source's  power,  its  interaction  in  any  known  environment  or  any  known 
surrounding  can  be  determined. 

Until  recently,  most  all  experimental  measurements  to  describe  an 
acoustic  source  were  accomplished  by  measuring  sound  pressure.   This 
quantity,  however,  was  the  result  of  a  source's  interaction  with  a 
specific  environment  and  utilizing  the  experimental  results  in  other 
applications  was  often  difficult.   One  technique  that  was  developed  to 
overcome  the  difficulty  in  relating  one  environment  to  other  was  to 
isolate  the  source  of  interest  and  to  simulate  a  free  field  environ- 
ment in  an  anechoic  chamber.   Here  the  radiated  sound  pressure  in  air 
was  measured,  then  related  directly  to  the  sound  intensity.   This 
technique  has  as  its  biggest  limitation,  the  availability  of  an  effec- 
tive anechoic  chamber,  particularly  for  use  in  low  frequency  analysis 
at  the  low  frequencies. 

Another  technique  is  to  measure  the  intensity  directly  utilizing  a 
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sound  intensity  probe.   The  theory  of  intensity  measurements  comes 
directly  from  two  fundamental  properties  of  acoustics.   Those  are  the 
conservation  of  momentum  and  the  definition  of  intensity.   Conserva- 
tion on  momentum  states  that: 

pdvT/dt   -  -  3p/3r  (2.1) 

which  can  be  solved  for  the  radial  component  of  the  particle  velo- 
city, vr 

vr  -  -  (1/p)  /  3p/3r  dt  (2.1a) 

A  finite  difference  approximation  is  used  to  calculate  the  pres- 
sure gradient  using  two  identical  microphones  by 


dp/dr   -  lim  (Ap/Ar)  *  (pA-pB)/Ar  (2.2) 

Ar-0 


where  A  and  B  correspond  to  channels  A  and  B  for  the  respective 
microphone  inputs  and  Ar  is  the  finite  microphone  spacing.  The 
instantaneous  pressure  is  taken  as  the  mean  value  of  the  two  pressure 
signals, 

P  *  (PA+PB>/2-  <2-2a> 

The  finite  difference  approximation  is  valid  provided  the  separa- 


.  14  - 

tion  distance,  Ar,  is  small  compared  to  the  acoustic  wavelength. 
Thus,  a  finite  difference  approximation  error  becomes  significant  at 
high  frequencies.  For  a  plane  sinusoidal  wave,  propagating  along  the 
axis  joining  the  two  microphones,  the  finite  difference  approximation 
assumes  the  free  field  phase  between  the  two  microphone  positions  to 
be  equal  to  kAr,  whereas  the  actual  measured  free  field  phase  is 
Sin(kAr).   Therefore,  the  measured  intensity,  I',  is  related  to  the 
actual  intensity,  I,  by 

I'/I  -  Sin(kAr)/kAr  (2.3) 

and  the  resulting  measured  intensity  underestimates  the  actual 
intensity  by  Le  at  high  frequency  by 

Le  -  10  log (Sin (kAr) /kAr)  .  (2.3a) 

All  the  intensity  measurements  in  this  thesis  were  obtained  with  a 
spacing  of  Ar  -  50mm.   Thus  the  finite  difference  error  for  this 
spacing,  at  the  upper  analysis  bandwidth  frequency  of  f  -  1600  Hz,  was 
Le  -  -  1.7  dB. 

From  the  definition  of  intensity  and  by  substituting  equations  2.2 
and  2.2a,  the  measured  sound  intensity  becomes 
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Ir  -  pvr  =  -  ((pA+pB)/2pAr]J(pA-pB)dt.  (2.3b) 

By  measuring  the  intensity  over  a  control  surface  which  surrounds 
the  source,  the  power  can  be  calculated  from: 

P  -  /  Ir  dA  (2.4) 

A 

The  advantages  of  measuring  intensity  directly  are:  (1)  the 

effects  from  diffuse  noise  contaminating  the  experimental  results  are 

eliminated  since  the  integral  of  this  noise  intensity  around  any 

closed  control  surface  is  zero,  (2)  the  sound  power  calculated 

describes  the  acoustic  source  In  a  form  which  is  readily  usable  to 

design  engineers,  (3)  measurements  can  be  made  in  the  source's  near 

field  and  (4)  there  are  no  restrictions  on  the  shape  and  size  of  the 

control  surface.   A  simple  yet  effective  analogy  to  describe  the 

rationale  in  measuring  sound  intensity,  thus  power  directly  follows: 

In  attempting  to  design  a  heating  system,  an  engineer 
solicits  specifications  from  two  companies.  One  com- 
pany describes  its  heating  system  by  stating  the 
electrical  power  and  the  thermal  efficiency  of  the 
unit.  The  other  company  describes  the  temperature  at 
various  locations  inside  an  experimental  enclosure  at 
steady  state  conditions. 

It  Is  obvious  in  this  example  that  describing  the  temperature 
field  within  one  environment  is  practically  useless  to  the  design 
engineer  whereas  knowing  the  power  and  thermal  efficiency,  the 
designer  can  determine  the  interaction  of  this  source  with  his  partic- 
ular environment. 
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In  the  experiments  which  follow,  every  attempt  was  made  to  uti- 
lize the  advantages  of  measuring  intensity  directly  and  simulating  a 
free  field  environment  by  taking  the  measurements  in  a  semi-anechoic 
chamber,  thus  minimizing  any  enclosure  interactions  which  could  affect 
the  radiated  sound  power. 

2.2   Ideal  Versus  Practical  Sound  Intensity  Analyzer 

2.2.1  Active  and  reactive  sound  fields 

Sound  fields  exist  in  two  parts:  an  active  part  which  transports 
sound  energy  by  virtue  of  the  fact  that  the  sound  pressure  and  the 
particle  velocity  are  in  phase  and  a  reactive  part  which  stores  sound 
energy  and  in  which  the  sound  pressure  and  particle  velocity  are  in 
quadrature.   A  sound  field  is  said  to  be  more  reactive  the  more  energy 
is  stored  relative  to  that  transported.   Sound  field  reactivity  is 
defined  as  the  difference  between  the  measured  sound  pressure  level, 
Lp,  and  the  maximum  intensity  level  L± . 

2.2.2  Ideal  sound  intensity  analyzer 

An  ideal  sound  intensity  analyzer  would  measure  only  the  active 
part  of  the  total  sound  field  and  would  indicate  sound  intensity  val- 
ues approaching  0  watts/m2   in  a  highly  reactive  (highly  diffuse) 
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sound  field. 

2.2.3   Practical  sound  intensity  analyzer 

A  practical  sound  intensity  analyzer  can  never  display  such  a  low 
intensity  level  due  to  the  small  but  ever  present  phase-mismatch 
between  the  two  microphone  channels.   It  is  this  phase  mismatch  which 
accounts  for  the  lower  limit  in  the  dynamic  range  of  a  practical  sound 
intensity  analyzer. 

2.3  Determining  the  Limits  of  Error  for  the  Measured  Intensity 

Errors  are  introduced  in  any  practical  intensity  measurement  by: 

(1)  the  influence  of  the  phase  mis -match  in  the  probe,  4>m   and 
the  analyzer,  <f>&, 

(2)  the  influence  of  the  phase  of  the  sound  field  of  known  reacti- 
vity, 4>f,    and 

(3)  using  a  particular  spacer  between  the  two  microphones. 

The  limits  of  these  errors  can  be  determined  using  a  phase- 
reactivity  nomogram  [see  reference  6].   This  nomogram  graphically 
determines,  <f>£   for  a  given  measured  reactivity  (L^-Lp)  and  a 
particular  spacer,  Ar,  according  to  the  relationship: 
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Lj_  -  (Lp+0.16)  -  10Log[1000/Ar  x  C/f  x  ^f/360]  (2.5) 

The  phase  of  the  sound  field,  <f>£,    together  with  the  phase  of  the  probe 
and  the  analyzer,  <j>£   •    <f>a    -    <f>m   and  <f>f   +  <f>&   +  ^m,  set  the  limits  of 
error  in  the  measured  phase  at  a  particular  frequency.   From  these 
phase  limits,  the  nomogram  provides  the  limits  of  error  in  dB  for  the 
measured  intensity. 

2.4  Sound  Intensity  Probe 

The  specific  details  on  the  operation  of  the  intensity  probe  used 
in  this  experiment  are  contained  In  reference  (7).   Briefly,  the 
intensity  probe  consists  to  two  B&K  phase  matched,  1/2"  condenser 
microphones,  type  4177/4165,  separated  by  a  spacer  50  mm  in  length  and 
on  a  common  axis.   One  microphone  is  a  reference  while  the  other 
microphone  provides  the  pressure  gradient  measurement. 

2.5  Intensity  Measurements 

2.5.1  Calculating  intensity 

Sound  intensity  measurements  were  analyzed  using  a  B&K  2032  Spec- 
trum Analyzer.   The  sound  intensity  was  calculated  according  to  the 
relationship: 
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I(k)  -  -  (1/u.Arp)  Im[GAB(k)] 


(2.6) 


where  Ar  is  the  microphone  spacing  in  meters,  A  corresponds  to 
Channel  A  which  is  the  reference  microphone  input,   B  corresponds  to 
Channel  B  which  is  the  second  microphone  input,  k  is  the  frequency 
index  and  p   is  the  density  of  air,  calculated  according  to  the  rela- 
tionship : 


p   -  0.348534  p/(273+T) 


(2.7) 


where  p  is  the  pressure  in  mbar  and  T  is  the  temperature  in  degrees 
centigrade . 


The  cross  spectrum,  G^(k)  ,  is  the  complex  ensemble  average  of  the 
complex  products  of  the  complex  conjugated  one-sided  instantaneous 
spectrum  GA(k)  and  the  one-sided  instantaneous  spectrum  Gg(k) . 


GAB<k>  "  GA(k)-GB(k)  (2.8) 

The  instantaneous  one-sided  spectrum  for  channel  A  or  B,  GA(k) , 
Gg(k)  was  calculated  according  to: 


GA(k)  - 


r  SA(k)      for  k  -  0 
2  SA(k)    for  1  <  k  <  N/2  -  1 
0         for  N/2  <  k  <  N-l 


(2.9) 
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where   N  is  the  number  of  samples  per  ensemble.   For  this  experiment, 
N  was  equal  to  1024.   SA(k)  is  the  two-sided  instantaneous  spectrum 
for  channel  A  which  is  determined  from  the  recorded  time  record,  a(n) 
according  to: 

Sa(k)  -?[w(n).a(n)]  (2.10) 

where  j    is  the  forward  discrete  fourier  transform  and  u>(n)  is  the 
time  record  weighting  function.   In  this  experiment,  a  Harming  window 
weighting  function  was  utilized. 


The  forward  discrete  fourier  transform,  j    ,  is  defined 


as : 


N-l 
SA(k)  -  1/N  2  a(n)«exp(-j  2wkn/N)  (2.11) 

n-0 


where  n  equals  the  time  index  times  the  sampling  interval,  nAt 
and  k  equals  the  frequency  index  times  the  frequency  resolution,  kAf 
The  frequency  resolution  for  this  experiment  was  2Hz  for  an  analysis 
bandwidth  of  1600  Hz. 

2.5.2  Measurement  techniques 

Two  techniques  can  be  used  to  measure  sound  intensity  using  an 
intensity  probe.   They  are: 
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(1)  continuous  spatial  averaging 

(2)  discrete  point  averaging 

Continuous  spatial  averaging  consists  of  sweeping  the  intensity 
probe  at  a  controlled  rate  over  each  surface  while  continuously  samp- 
ling and  averaging  data.   This  technique  produces  the  most  accurate 
spatially  averaged  results  but  suffers  slightly  in  repeatability. 
Several  repeated  runs,  often  four  or  more,  were  required  to  obtain  the 
required  repeatability  of  total  sound  power  (within  0.5  dB  on  at  least 
three  runs).   Figures  2.1,  2.2  and  2.3  are  the  results  of  three  such 
repeatability  experiments  with  0.5  dB  or  less  variation  in  total  sound 
power.   Shown  is  the  sound  power  spectral  level  from  a  rectangular 
orifice,  4  1/2"  by  1  2/3",  and  an  analysis  bandwidth  of  2  Hz. 
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Figure  2.1:  Louver  11  sound  power  spectra  repeatability.  Shown  are  four  succes- 
sive spectral  measurements  of  sound  power. 
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Fiflure  2.2:  Orifice  01  sound  power  spectra  repeatability.  Shown  are  three  suc- 
cessive spectral  measurements  of  sound  power. 
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Figure  2.3  Orifice  01  sound  power  spectra  repeatability.  Shown  are  five  succes- 
sive spectral  measurements  of  sound  power. 
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Discrete  point  averaging  consists  of  dividing  your  control  surface 
Into  numerous  discrete  areas  and  measuring  the  sound  Intensity  at  the 
center  of  each  discrete  area  (point)  and  summing  over  all  areas.   This 
technique  offers  maximum  repeatability  but  at  reduced  spatially- 
averaged  accuracy.  The  two  techniques  where  both  utilized  initially 
and  compared.   The  repeatability  of  the  total  radiated  power  using  the 
discrete  point  technique,  with  45  points,  was  0.3  -  0.5  dB  between 
successive  runs.   These  repeatability  values  were  selected  as  the 
criteria  for  the  continuous  spatial  averaging  technique.   The  results 
showed  that  the  two  techniques  gave  similar  sound  power  spectrums  and 
similar  total  sound  power  levels  within  an  accuracy  of  0.3  dB.   For 
this  reason,  the  continuous  spatial  averaging  technique  was  selected 
due  to  its  increased  accuracy  and  reduced  sampling  time. 

The  control  surface  was  constructed  of  1/4"  balsa  wood  enclosing  a 
total  surface  area  of  0.542  square  meters.   Intensity  spectra  over  the 
frequency  range  of  52-1600  Hz  were  taken  and  recorded  on  a  B&K  2032 
Spectrum  Analyzer.   Microphone  calibration  was  performed  using  a  pis- 
tonphone  and  daily  temperature  and  barometer  measurements  were  taken 
to  calculate  the  air  density.   To  ensure  that  the  two  microphone  were 
indeed  phase  matched,  their  phase  was  measured  in  a  random  noise  field 
and  compared.  The  microphones  were  placed  at  the  open  end  of  a  long 
cylinder,  equal  distances  from  a  speaker  placed  at  the  opposite  end  of 
the  cylinder.   This  apparatus  was  similar  in  design  to  a  calibration 
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pistonphone,    type   4200.    Figure   2.4   shows    the   phase   difference   measured 
over   a   frequency   range   of   2560  Hz   and  clearly   demonstrates    their   com- 
patibility over   the    frequency  band  of   interest,    52    -    1600  Hz. 
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Figure  2.4:     Phase  difference  between  two  microphones  and  their  pre-anplif iers 
(channels  A  and  B) 


In   Section  5.2,    the   errors    in  dB   resulting   from   these   phase  mis- 
matches  are   calculated  and   shown   to  be   on   the    same   order   of  magnitude 
as    the    inherent  statistical   fluctuations. 
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2.5.3   Power  spectrum  contamination 

Data  below  52  Hz  were  ignored  due  to  contamination  by  the  1,0,0 
room  mode  (38  Hz)  of  the  semi-anechoic  chamber.   Also  at  this  low 
frequency,  the  effectiveness  of  the  anechoic  treatment  was  reduced. 
The  peak  at  462  Hz  and  its  harmonic  at  924  Hz  which  is  evident  in  all 
the  spectra  was  due  to  the  resonance  of  the  cross -duct  mode  of  the 
test  section  and  not  at  all  associated  with  the  orifice  source.   For- 
tunately, the  sound  level  and  narrow  bandwidth  associated  with  this 
modal  resonance  contributes  insignificantly  to  the  total  sound  power 
caused  by  the  orifice  source.  Frequencies  above  1600  Hz  contributed 
nothing  to  the  total  sound  power  and  the  spectral  shape  followed  the 
same  trend  as  from  1200  -  1600  Hz. 
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3 .   THEORY 
3.1  Introduction 

The  development  of  a  theoretical  model  to  determine  the  scaling 
laws  for  orifice/louver  flow  noise  was  first  accomplished  in  1972  by 
Ffowcs  Williams  [8].   He  treated  the  model  problem  of  flow  noise  gen- 
erated by  a  body  of  turbulence  near  an  Infinite,  rigid,  thin  screen 
with  circular  orifices.  His  assumptions  were:  a)  the  size  of  the  ori- 
fice is  much  smaller  than  the  distance  between  them,  b)  variations  in 
pressure  within  individual  apertures  are  neglected.   While  these 
assumptions  are  rather  restrictive  in  their  application  to  many  indus- 
trial louvers,  the  theory  does  provide  a  fundamental  result:  that  the 
sound  radiated  by  orifices  and  louvers  can  be  directly  related  to  the 
mass  flow  driven  through  the  aperture  by  the  aerodynamic  field. 

This  fundamental  result  was  used  by  Nelson  [9]  to  develop  a  model 
for  flow  noise  generated  over  practical  perforated  screens  with  circu- 
lar orifices.   Due  to  the  normal  component  of  the  turbulent  fluctua- 
tions above  the  louver,  flow  is  driven  through  the  orifice  resulting 
in  a  fluctuating  volume  flow  through  the  upper  surface.   This  fluctu- 
ating volume  flow  constitutes  an  equivalent  monopole  source  at  the 
upper  surface.   Since  the  volume  flow  into  the  upper  surface  of  the 
orifice  creates  an  equivalent  outflow  through  the  lower  surface,  there 
is  an  equal  and  opposite  equivalent  monopole  at  the  lower  surface. 


-  27  - 

The  net  field  which  radiates  out  of  the  louver  and  into  the  closed 
acoustic  space  is  thus  the  combination  of  the  lower  monopole  source 
and  the  contribution  from  the  upper  monopole  source  which  is  dependent 
on  the  degree  of  acoustical  transparency  of  the  louver.   If  the  louver 
is  acoustically  opaque,  then  only  the  lower  monopole  source  will  radi- 
ate into  the  acoustic  space.   If  the  louver  is  acoustically  transpar- 
ent, then  the  upper  monopole  source  will  radiate  into  the  space,  hav- 
ing a  cancelling  effect  on  the  field  of  the  lower  monopole  source. 
This  would  result  in  a  net  radiated  field  which  exhibits  properties  of 
a  dipole  source. 

Nelson's  equivalent  source  model  gave  an  excellent  description  of 
the  flow  noise  produced  by  perforated  facing  materials  used  in  prac- 
tice even  though  he  assumed  a  rather  simplistic  form  for  the  driving 
aerodynamic  velocity.  His  assumptions  were:  a)  the  wavelength  for  the 
sound  is  much  larger  than  the  orifice  radius,  b)  there  are  many  aper- 
tures per  acoustic  wavelength  and  no  interactions  of  the  flow  occurred 
between  adjacent  apertures  which  can  modify  the  aperture's  conducti- 
vity.  In  the  sections  that  follow,  the  theory  developed  by  Nelson  and 
Ffowcs  Williams  will  be  utilized  to  model  the  flow  noise  created  by 
turbulent  flow  over  orifices  and  louvers.  Once  the  equivalent  source 
model  has  been  developed,  it  will  be  used  to  collapse  the  experimental 
sound  power  spectra  thus  providing  a  verification  to  the  theory's  val- 
idity.  The  coordinated  system  used  for  this  analysis  is  shown  in 
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Figure  3.1 


Figure  3.1:  Coordinate  system  for  analysis 
3.2  The  Equivalent  Source  Model  for  Orifices  and  Louvers 

Ffowcs  Williams  [8]  and  Leppington  and  Levine  [10]  developed  an 
expression  for  the  reflection  coefficient,   R,  for  a  plane  wave  of 
single  frequency,  u,  incident  at  a  polar  angle,  $,    on  an  infinite 
louver.   For  low  frequencies  and  orifice  spacings  that  are  large  com- 
pared to  the  orifice  size,  the  reflection  coefficient,  R,  is 


R  -  l/(l+4aiNAcos0) 


(3.1) 


where  a   is  the  orifice  radius,   N  the  number  of  orifices  per 
unit  area  of  louver  and  k   is  the  acoustic  wavenumber,  k  -  u>/Co . 

The  field  produced  by  the  lower  monopole  at  y^  can  be  determined 
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using  the  principle  of  reciprocity.   The  field  at  y^  due  to  a  distant 
source  at  x^,  is  given  by  p^(l+R)  where  p^  is  the  incident  plane  wave 
striking  the  lower  surface  of  the  louver.    Thus  one  can  write 

P(yic>  "  <1+R>  Q<xk>  exp(ik|yk-xk|  )/4^r|yk-xkl  (3-2) 

where  Q  is  the  strength  of  the  distance  monopole  at  x^.   By 
reciprocity,  the  field  at  x^  due  to  the  monopole  at  y^  is  given  by 
interchanging  the  source  and  receiver  position  vector  in  equation  3.2. 
Equation  3.2  can  be  generalized  to  related  the  spectral  density  of  the 
acoustic  pressure  field,  SpCx^.w)  to  the  spectral  density  of  the 
source  SqCy^.w)   such  that 


2.2 


Sp(xk,«)  -  |l+R|z.SQ(Yk,w)/16^r 


(3.3) 


where  r*  -  Ix^-y^p.   Figure  3.2  is  a  plot  of   |Rp  against  the  rati< 
(4aN/kcos0). 


i*i*  •<#».#'«•*  i* 


K.*1imH 


Figure  3.2:  The  infinite  screen  reflection  coefficient  as  a  function  of  the  ratio 
(4aN/kcos0>.  From  [9,  Nelson,  P. A.,  1982]. 
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From  this  plot,  it  can  be  seen  that, 

|R|2  -  1  for  (4aNAcos0)  s   0.3. 

This  occurs  for  open  area  ratios  on  the  order  of  0.01  (see  refer- 
ence 11).  For  all  the  orifice  and  louvers  tested  in  this  study,  the 
open  area  ratio  is  less  than  0.005  so  that  |r|2  -  1.   To  verify  this 
result  for  the  reflection  coefficient,  since  the  assumption  of  many 
apertures  per  acoustic  wavelength  is  not  valid  for  the  single  orifice, 
a  separate  analysis  was  performed.   Here,  in  two  dimensions,  the  plane 
wave  of  single  frequency,  w,  is  incident  normal  on  an  infinite,  rigid 
surface  with  a  single  opening,  of  dimension,  2b.  Again  the  principle 
of  reciprocity  is  used  to  predict  the  radiated  pressure  from  the 
single  orifice.   Using  matched  asymptotic  expansions,   for  b/A  «  1, 
the  lowest  order  solution  for  the  reflection  coefficient  is  again, 
|r|2  -  1.  The  details  of  this  analysis  are  shown  in  Appendix  1. 


Thus, 


.2^2 


Sp(xk,to)  -  SQ(yk,c)/4*^  (3.4) 


The    instantaneous   source   strength  q(yk,t)    can  be  written  in  terms 
of  the  weighted   integral  of  the  unsteady  normal   component  of  velocity 
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over  the  area  of  the  aperture.   To  take  advantage  of  the  homogeneous 
and  stationary  properties  of  the  velocity  in  the  plane  of  the  wall,  we 
can  write 

q(yk»c>  -  -  i(^J  /v(yk.t)K<sk-yk>dA(sk)  <3-5) 

CO  00 

where  yk  is  the  vector  position  of  some  reference  point  on  the  aper- 
ture surface  and  K(sk-yk)  is  the  weighting  function,  i.e.,  the 
function  which  converts  the  integral  over  the  aperture  area,  sk  - 
sl»s2»  to  an  integral  over  all  space,  defined  as 

K(sk-yk)  -  1  for  sk-yk  inside  of  the  aperture  area,  A 

-  0  for  sk  -  yk  oustide  A 
Following  Corcos  [12]  the  spectral  density  of  the  source  strength, 
Sn(yk,w)  can  be  defined  by 

SQ(yk.")  "  P2»2S   JSvv(2k-sk,w)K(zk)K(sk)dA(zk)dA(sk)        (3.6) 

00  00 

where   zk   is  a  dummy  space  vector.   Define  ck  by 

€k  "  zk  •  sk 

thus 

sQ(yk'w)  "  P1"1  S   /Svv(€k,w)K(sk)K(sk+ek)dA(sk)dA(ek) 

co  oo 

-  P2o>2   J   Svv(ck,u>)tf(ck)dA(ek)  (3.6b) 

00 

where  the  function,   0(ek),  defined  as  the  orifice  response  function, 


-  32  - 

is  determined  wholly  by  the  geometry  of  the  aperture  and  can  be  evalu- 
ated as 

00 

*(<k)  "  jK(sk)K(sk+€k)dA(sk)  -  1  x  area  of  overlap 


-  (L-|e3|)<b-|e1|)  for  |e3|  <  L,  |ei|  <  b 

-  0  elsewhere 

for  the  rectangular  shaped  orifice  shown  below 


(3.6c) 


/ 

\ 

s 

< 

Finally,  we  assume  the  cross-spectral  density,  Svv(ek>w)»  to  have 
the  product  form 


iK«, 


Sw(<k,u>)  "  Sv(o>)A(K€1)B(Kc3)e1N£l 


(3.7) 


where  Sv(u>)  is  the  spectral  density  of  the  normal  velocity  fluctua- 
tions,  K  is  the  convective  wavenumber  given  by  K  -  <*>/Uc  where  Uc  is 
the  convection  speed  of  the  boundary  layer  eddies  and  A(Kc^)  and 
B(K«3)  are  the  amplitude  functions  for  the  normalized  longitudinal  and 
lateral  cross-spectra,  respectively.  This  form  is  similar  to  the 
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expression  given  by  Corcos  [12]  for  the  cross  spectral  density  of  the 
boundary  layer  pressure  fluctuations.  The  longitudinal  cross-spectra 
of  the  eddies  as  they  convect  in  the   y^  direction  also  contains  the 
phase  term  expCiKc^).   ACKt^)  is  plotted  in  Figure  3.3  as  a  function 
of  its  argument  for  various  ratios  of  separation  distance  to  momentum 
thickness  (r^/fi  )  for  wall  pressure.  Since  the  wall  pressure  is 
related  to  the  volume  integral  of  the  normal  component  of  the  velo- 
city, wall  pressure  cross-spectra  data  can  be  used  to  predict  the  lon- 
gitudinal and  the  lateral  cross -spectra  of  the  normal  component  of 
particle  velocity.   Comparing  the  wall  pressure  correlation  length 
scales  from  [12]  to  the  velocity  correlation  length  scales  from  [17], 
we  find  that  both  the  longitudinal  and  lateral  length  scales  agree  to 
within  an  order  of  magnitude,  thus  supporting  the  previous  statement. 


Figure  3.3:  Amplitude  function  for  longitudinal  cross- spectral  density  of  wall 
pressure.  From  [16,  B8N  report  No.  4110366,  1966] 
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For  the  various  orifice  and  louver  samples  tested,  at  a  maximum 
value  for  r^/6   equal  to  1,  A(K«^)  varies  from  a  maximum  value  of  1  to 
a  minimum  value  of  0.6  for  orifice  02  at  U-20m/s  and  f-250  Hz  (the 
frequency  at  which  the  spectral  level  is  lOdB  down  from  the  peak 
level).   Throughout  most  of  the  samples  tested,  A(Kf^)  -  1  over  the 
frequency  band  of  interest.   Thus  the  longitudinal  correlation  lengths 
were  typically  several  times  the  length  of  the  orifice  aperture.   Thus 
the  assumption  is  made  that  no  decay  in  the  eddy  structure  occurs  in 
the  longitudinal  (y^)  direction. 

The  normalized  lateral  cross-spectra  in  the  y$   direction  is  given 
by  the  term  B(K«3)  which  is  assumed  unity  for  no  separation  and  to 
decay  very  rapidly  with  increased  separation.   To  justify  this  assump- 
tion, B(K«3)  is  plotted  in  Figure  3.4  as  a  function  of  its  argument 
for  various  ratios  of  x-^/S    .   From  this  figure,  the  typical  correla- 
tion length,  i,  is  calculated  to  be  on  the  order  of  2  cm.   The  ratio 
of  the  lateral  correlation  length,  i,  to  the  aperture  length,  L,  is 
0.18  thus  supporting  the  initial  assumption  of  rapid  lateral  decay. 
The  correlation  lengths  are  determined  at  the  point  where  the  ampli- 
tude function  equals  0.4  and  at  the  frequency  corresponding  to  the 
peak  power  spectral  level. 
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Figure  3.4:  Amplitude  function  for  lateral  cross -spectra I  density  of  wall  pressure 
From   [16,    BBN  Report  No.    4110366,    1966] 


Substituting  equation  3.7  into  equation  3.6b  and  evaluating  the 
integral  over  the  rectangular  shaped  orifice,  assuming  SyCy^.w)  to  be 
uniform  over  the  orifice  area  gives: 


SQ(yk.")  "  P 2»2   Sv(yk,u>).Li.(2/K2)(l-cos(Kb)) 


(3.8) 


where  L  is  the  orifice  length  in  the  y$   direction.   See  Appendix  (3) 
for  the  details  of  the  integration  for  the  rectangular  shaped  orifices 
and  reference  [9]  for  the  results  for  circular  shaped  orifices. 
Equation  3.8  states  that  the  monopole  spectral  density  is  directly 
related  to  the  normal  velocity  spectral  density. 
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3.3  Scaling  Laws  for  the  Radiated  Sound  Power 

Having  determined  the  form  for  the  source  strength,  we  can  now 
develop  the  scaling  laws  for  the  radiated  sound  power  from  various 
rectangular  shaped  configurations.   Assuming  that  each  orifice  radi- 
ates sound  incoherently,  the  total  source  strength  for  N  apertures 
is  simply  N   times  the  individual  source  strength  for  each  aperture 
defined  by  equation  3.8.   Substituting  equation  3.8  into  equation  3.4 
gives : 

2 
Sp(xk,o)  -  Np2U.Sv(yk,w)(l-cos(Kb))Uc/2*2r2  (3.9) 

The  assumption  is  made  that  the  convection  velocity,  Uc,  is  equal 
to  one  half  the  frees tream  velocity,  Uw.   Define  the  turbulent 
intensity,  I,  and  the  Strouhal  number,  St(i)  and  St(b)  by: 

(I2)^  -  1/U.  •  /  Sv(yk,u>)  du>  (3.10) 

St(i)  -  fi/U..   St(b)  -  fb/U.  (3.11) 

Finally,    integrating  equation  3.9  over  the  bandwidth  Ao>  with   the   above 
expressions   substituted,    yields 


5 

2n.     ,1     /Of,r2 


(pMxk,t))Aw  -     Np^L.(l-cos(Kb))St(i)(I^)AwU<D/2fxrz  (3.12) 
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For  spherical  waves,  the  sound  power,  P,  can  be  determined  by: 


(P(t))Aw  -  NpL  (l-cosCKb^StUXI2)^  U./fCo  (3.13) 

Equation  3.13  is  the  final  result  for  the  flow  noise  created  by  tur- 
bulent flow  over  orifices  and  louvers  with  rectangular  shape  as 
described  in  Section  (4.2).   The  controlling  parameters  for  scaling 
sound  power  are  thus: 

(1)  Strouhal  number  based  on  the  lateral  correlation  length,  i 
and  aperture  transverse  dimension,  b 

(2)  the  normal  component  of  the  Turbulent  Intensity 

(3)  the  fifth  power  of  the  freestream  velocity 

(4)  the  shape  factor  (l-cos(Kb)) 

In  Section  (5.10),  this  model  is  used  to  collapse  the  experimental 
sound  power  data  for  orifices  02  and  03  and  louver  LI.   It  is  clear 
from  the  collapse  that  equation  3.13  does  Indeed  identify  the  predom- 
inant controlling  parameters  for  rectangular  shaped  orifice  and  louver 
flow  noise.   Note,  for  circular  shaped  apertures,  the  above  scaling 
factors  remain  the  same  except  the  shape  factor  is  replaced  by  the 
function,  H^(Kd)  where  H^  is  a  Struve  function  of  the  first  kind  and 
d   is  the  aperture  diameter.   See  reference  9  for  the  details  in  cal- 
culating equation  3.8  for  a  circular  aperture. 
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4.   EXPERIMENTAL  APPARATUS 

4.1  Seml-anechoic  Wind  Tunnel 

The  experiments  were  carried  out  using  the  low  turbulence,  sub- 
sonic wind  tunnel  of  the  Acoustics  and  Vibration  Laboratory  at  MIT. 
The  wind  tunnel  details  are  shown  in  Figure  (4.1)  with  the  basic  con- 
struction characteristics  described  by  Hanson  [13]  and  modifications 
described  by  Shapiro  [14] .   The  wind  tunnel  is  of  open  circuit  design 
with  the  test  section  enclosed  in  a  blockhouse.   The  tunnel  consists 
of  a  set  of  flow  straighteners,  screens,  a  settling  chamber  and  a  20:1 
contraction  leading  to  a  square  test  section  38cm  by  38cm.   The  block- 
house walls,  floor  and  ceiling  were  treated  with  urethane  foam  to  pro- 
vide a  semi-anechoic  chamber  within  which  acoustic  intensity  measure- 
ments were  taken.   Figure  (4.2)  details  a  cut-away  front  and  top  view 
of  the  chamber  and  test  section.   The  acoustic  treatment  consisted  of 
a  base  of  10cm  foam,  then  a  set  of  10cm  foam  squares  spread  randomly 
and  finally  a  set  of  2cm  foam  squares  also  spread  randomly.   The 
effectiveness  of  the  coating  was  determined  by  [14]  using  a  horn 
driver  and  measuring  the  sound  field  decay  with  distance.   The  ane- 
choic  treatment  was  effective  between  200  Hz  and  2000  Hz  and  only  par- 
tially effective  at  100  Hz. 
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Figure  4.1:  Wind  Tunnel  Facility 
4.2  Orlf ice/Louver  Samples 

Nine  separate  orifice/louver  samples  were  tested.   All  samples 
were  made  of  plexiglass  and  screwed  into  the  sample  mounting  window 
which  was  mounted  flush  into  the  side  to  the  test  section  as  shown  in 
Figure  (4.2).   The  separate  samples  tested  were: 


(1)  0-1:   a  4  1/2"  (height)  by  1  2/3"  (width)  rectangular  ori- 

fice, 3/4"  thick  with  straight  sided  walls 

(2)  0-2:   the  above  made  of  1/4"  thick  plexiglass 
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Figure  4.2:  Experimental  set-up  within  the  wind  tunnel,  front  and  top  views 
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(3)  0-3:   a  4  7/16"  by  27/64"  rectangular  orifice,  1/4"  thick 

with  straight  sided  walls 

(4)  0-2A:  0-2  with  40  degree  slanted  walls,  slanted  away  from  the 

flow  to  reduce  the  downstream  edge  impact  angle 

(5)  0-3A:  0-3  with  40  degree  slanted  walls 

(6)  L-l:   4  rectangular  orifices,  each  4  7/16"  by  27/64",  separ- 

ated by  1/4",  1/4"  thick 

(7)  L-2:   54  circular  orifices,  each  27/64"  in  diameter,  separ- 

ated by  0.6328"  between  centers,  1/4"  thick 

(8)  L-1A:  L-l  with  40  degree  slanted  walls 

(9)  L-2A:  L-2  with  40  degree  slanted  walls. 

The  louver  configurations  were  designed  to  maintain  the  same  total 
open  area  as  the  single  orifices,  0-1  and  0-2.  The  40  degree  slanted 
walls  were  slanted  to  provide  a  blunt  trailing  edge  and  subsequently  a 
sharp  leading  edge . 

4.3  Test  Section 

The  test  section  was  made  from  3/4"  finished  plywood,  38cm  by  38cm 
square  in  cross  section  and  was  located  on  the  centerline  of  the  cham- 
ber.  The  test  section  was  sealed  and  the  ends  were  further  covered  by 
2"  of  ure thane  foam.   The  background  power  levels,  measured  with  the 
orifice  covered,  were  sufficiently  low  compared  to  the  power  levels 
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with  the  orifice  exposed  (see  Section  5.3)  to  ensure  signal  to  noise 
ratios  of  at  least  10  dB  over  the  frequency  range  of  interest  and  at 
all  operating  speeds. 

4.4  Turbulent  Intensity  Measurements 

Turbulent  Intensity  is  defined  by  the  relationship: 

TI  -  20  log  (/uVu.)  (4.1) 

Turbulent  fluctuations  were  measured  using  a  Dantec  Constant  Tempera- 
ture Anemometer  with  a  5  micron  diameter  tungsten  wire.   The  output  of 
the  anemometer,  being  the  total  voltage  E,  was  inputted  directly  to 
the  B&K  2032  Spectrum  Analyzer.   Here  the  AC-coupled  autospectrum, 
Ge(u) ,  was  measured.   The  mean-squared  fluctuating  velocity  autospec- 
trum, Gu(w)  was  determined  from: 

Gu(w)  -  (|H|2)  GE(w)  (4.2) 

where  H  is  the  transfer  function.  See  Appendix  2  for  the  details  in 
determining  H.  Turbulent  Intensity  spectra  were  measured  at  all  four 
flow  speeds  for  orifice  0-2  and  louvers  L-l  and  L-lA.  The  position  of 
the  U  probe  was : 
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(1)  two-tenths  of  an  inch  into  the  boundary  layer  directly  above 
the  orifice  at  both  the  upstream  and  downstream  edge  (for  the 
louver,  the  streamwise  measurements  were  taken  in  the  center 
of  the  furthermost  upstream  and  downstream  orifice,  and 

(2)  half-way  into  the  throat  of  the  orifice  and  at  the  same 
streamwise  location  as  described  above. 


-  uu 


5.   EXPERIMENTAL  RESULTS 
5.1  Introduction 

Figure  5.1  shows  the  geometry  of  the  problem.   Air  flows  at  free 
speed  Ug,  over  an  orifice  or  louver  which  has  been  inserted  flush  into 
a  long  flat  walled  duct.   The  orifice  was  either  rectangular,  with 
streamwise  dimension  b  and  length  L  or  circular  with  diameter  d. 
The  acoustic  radiated  power  from  the  orifice/louver  was  measured  at 
free  stream  velocities  of  20,  30,  40  and  50  meters  per  second.   The 
acoustic  power  was  determined  by  sweeping  the  intensity  probe  over  a 
control  surface  (a  five  sided  boxed  frame  surrounding  the  orifice) 
then  multiplying  the  resulting  spatially  and  temporally  averaged 
intensity  by  the  total  surface  area  of  the  control  surface.   Five  ori- 
fices and  four  louvers  were  tested  as  described  in  Section  4.2. 

u 


et«c 


^ mmnnr    nninnnp^ 


32.9^ 


\ Cot\\roL 

i^f  jet 

«3ch  •»!<!*  dimension  3£.9cm 
Figure  5.1:  Schematic  of  orifice  and  probe  geometry.  Air  flows  at  freestream  velo- 
city, U,  over  an  aperture  of  streamwise  dimension,  b.  Sound  intensity 
is  measured  by  probe  (microphone  spacing,  Ar)  while  continuously  sweep- 
ing control  surface. 


45 


To  determine  the  contribution  made  solely  by  the  orifice/louver 
source,  two  additional  experiments  were  conducted.   One  experiment 
consisted  of  measuring  the  acoustic  power  radiated  by  the  duct  when 
the  orifice  was  covered.   The  second  experiment  consisted  of  measuring 
the  power  flow  through  the  orifice,  in  the  absence  of  free  stream 
flow,  due  to  an  artificially  created  acoustic  field  in  the  duct.   This 
field  was  created  using  a  speaker  in  the  duct  and  was  equal  in  magni- 
tude to  that  caused  by  the  turbulent  boundary  layer  at  U  -  50  m/s . 
This  experiment  determined  the  "open-window"  contribution  to  the 
overall  sound  power.   Finally,  directivity  measurements  were  taken  for 
orifice  0-2  and  louver  L-2A  at  U  -  25  m/s  and  the  respective  polar 
plots  produced. 

5.2  Reactivity  of  the  Sound  Field  and  Sound  Intensity  Error  Limits 

The  reactivity  of  the  sound  field  in  the  blockhouse  (the  closed 
acoustic  space)  was  measured  for  various  orifice  and  louver  samples 
(see  Figures  5.2  and  5.3).   The  Table  below  lists  the  measured  dif- 
ference between  the  total  sound  pressure  level,  L-  and  the  total  sound 
intensity  level,  L^,  for  the  orifice  02A  and  the  louver  L2A  at  a 
radial  distance  of  2  feet. 
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Sample 


02A 


L2A 


j         Reactivity  <Lpl,> 

U„(m/s)                 | 

0.8 

30 

1.0 

50 

0.9 

30                     | 

1.0 

50 

Table  5.1:  Measured  Sound  Field  Reactivity 


»3;  -  - 
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Figure  5.2:  Sound  field  reactivity  for  orifice  02A.  Shown  are  the  sound  inten- 
sity and  sound  pressure  spectral  levels.  The  difference  between  the 
two  is  defined  as  the  sound  field  reactivity. 
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Frequency    (Hz) 


Figure  5.3:     Louver  L2A  sound  intensity  spectral   levels  versus  sound  pressure 
spectral   levels  (Lj  vs.  I  ).     The  difference  between  L|  and  L_  is 
defined  as  the  sound  field  reactivity. 


From  Figure   2.4,    the  measured  phase -mismatch  for   the  sound   inten- 
sity probe   is: 


Frequency  (Hz) 

63 

500 
1000 
2000 


4m  (degrees) 

0.8 
0.2 
0.5 
0.9 


The  maximum  phase -mismatch  between  channels  A  and  B  for   the  Ana- 
lyzer  is   0.3   degrees.      With  a  spacer  value   of  50mm  and  the   above   sound 
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field  reactivity,  the  phase  reactivity  nomogram  of  reference  6  gives 
the  following  value  for  the  maximum  phase  of  the  sound  field,  between 
the  two  microphone  positions. 

Frequency  (Hz)  *f  (decreet) 

63  2 

500  16 

1000  25 

2000  60 

Therefore  the  limits  of  the  phase -mismatch  are: 


Frequency  (Hz)  *f*a"*B  threes)  *f4*a**«  (deQrees) 

63  0.9  3.1 

500  15.5  16.5 

1000  24.2  25.8 

2000  58.8  61.2 


The  above  total  phase -mismatch  values  are  the  limits  of  error  in 
the  measured  phase  at  the  frequency  listed.   The  corresponding  limits 
in  dB  for  the  measured  intensity  from  the  nomogram  are: 

Frequency  (Hz)  Error  limits  (dB) 

63  ♦  2.5 

500  ♦  0.5 

1000  ♦  0.25 

2000  ♦  0.125 


The  corresponding  statistical  standard  deviation,  c,  of  the  mea- 
sured spectral  levels  about  the  mean  value  is: 
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<  -  4.34  /  7BTa  ,  (dB)  (5.1) 

where   B   is  the  filter  bandwidth  and  Ta  the  total  time  record  length. 
Substituting  the  appropriate  values  for  B(2  Hz)  and  Ta  (150  sec) 
gives  e   -   0.25  dB.   Thus  the  limits  of  error  in  measuring  the  sound 
intensity  are  approximately  equal  to  the  inherent  statistical  fluctua- 
tions (for  frequencies  above  100  Hz). 

5.3  Length  and  Frequency  Scaling 

The  three  significant  length  scales  for  this  experiment  were: 

(1)  the  streamwise  orifice  dimension,  b 

(2)  the  turbulent  boundary  layer  thickness,  S,    above  the  orifice, 
and 

(3)  the  lateral  correlation  length,  i. 

Table  5.2  provides  the  applicable  values  for  b,  6   and  i  used  in 
this  experiment. 

The  significant  frequency  scales  were  the  non-dimensional  strouhal 
numbers  St(b)  and  St(i)  defined  by: 


50 


St(b/2)   -  fcb/U*, 

St(i)        -  wi/Uc   -  4*fcf/U0 


(5.2) 


|  Orifice/ 

f/b 

t/b             | 

Louver 

b(cm) 

U(m/s) 

I  (cm) 

<5(on) 

for  0-2 

for  0-3 

0-1 

4.23 

20 

2.4 

4.19 

0.99 

3.95      j 

0-2,/ 

4.23 

30 

1.8 

3.66 

0.91 

3.64 

|    0-2A 

0-3,/ 

1.07 

40 

0.9 

3.65 

0.86 

3.44 

|    0-3A 

L-1,/ 

1.07 

50 

1.4 

3.49 

0.83 

3.29 

|    L-1A 

1-2./ 

1.07 

50 

2.1 

3.49 

|    L-2A 

Full 

2.50 

16.5 

4.70 

1.88 

Scale 

Note:  Full  scale  assumes  a  louver  positioned  15ft  from  the  leading  edge  of  a 

flat  surface  in  seawater  at  58  F.  The  louver  consists  of  aultiple  rectan- 
gular shaped  openings,  each  of  transverse  dinension,  b=2.5an.  From  Figure 
5.60  at  S/b  *   1.88,  the  Strouhal  #  «  0.19.  Therefore  the  full  scale 
frequency,  fc,  is: 

fc  «  (0.19X16.5)/  (0.025)  «  125  Hz 

Table  5.2  Experimental  Length  Scales 


where  fc  is  the  center  frequency  of  the  1/3  octave  band  with  the  peak 
sound  power  level.  Figure  5.4  is  a  graph  of  fc  versus  U^  for  orifice 
0-1  at  a  constant  b.   The  data  points  reflect  several  separate 
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experiments  measuring  the  sound  pressure  level  at  radial  distances  of 
9,  18,  24  and  36  inches  from  the  orifice.   Figure  5.5  plots  fc  versus 
Uw  for  orifices  0-1,  0-2,  0-3  and  louver  L-l  and  L-2.   Clearly  there 
exists  a  linear  relationship  (or  nearly  linear)  between  fc  and  U^  for 
a  particular  transverse  dimension  b. 
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Figure  5.4:     Frequency/velocity  relationship  for  orifice  01.     Shown  are  four  sep- 
arate experiments  for  various  radial  distances,  r. 


Table   5.3   lists   the  non-dimensional   Strouhal  numbers   for   the  above 
orifices/louver  as  well  as   the  non-dimensional   length  scale,      S/b   for 
U  -  20m/s.        The   strouhal  number,    St-fb/U,    has  been  defined  to  equal 
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that  recommended  by  Demetz  and  Farabee  [5]  to  compare  results.   Figure 
5.6  is  the  graph  of  the  data  in  Table  5.3  along  with  the  results  of 
Demetz  and  Farabee  [5]  for  TBL  flow  over  orifices  enclosed  by  solid 
cavities . 
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Figure  5.5:     Frequency/velocity  relationships  for  orifices  01,  02,  03  and  louvers 
L1.L2 
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|  Orifice/ 

St(b)  ■  fb/U 

Moo- dimensional 

Louver 

Strouhal  # 

length  scale  i/b 

|     0-1 

0.17 

1.22 

0-2 

0.23 

1.22 

0-3 

0.12 

4.79 

L-1 

0.076 

4.79 

I    L*2 

0.048 

4.79 

Table  5.3:  Non-dimensional  Frequency  and  length  Scales  at  U  =  20  m/s. 
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Figure  5.6:     Strouhal  number  versus  boundary  layer/aperture  width  ratio 
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The  decreasing  trend  in  the  Strouhal  number  with  increasing  6/b 
for  the  two  cases  agrees  reasonably  well.   Having  a  knowledge  of  the 
boundary  layer  thickness,  one  can  enter  Figure  5 . 6  at  the  appropriate 
value  for  5/b  and  determine  the  Strouhal  number  for  orifice  or  louver 
flow  noise.   From  the  Strouhal  number  and  the  geometry  of  the  apera- 
ture(s),  one  can  predict  fc. 

5.4  Structural  Radiation  Due  to  TBL  Excitation 

To  determine  the  contribution  from  structural  radiation  caused  by 
the  excitation  from  the  TBL,  an  experiment  was  conducted  at  all  free 
stream  velocities  with  the  orifice/louver  covered.   Identical  sound 
power  spectrums  were  recorded  as  described  in  Section  4.1.   The 
resulting  background  power  spectrums  levels,  Figures  5.7,  5.8,  5.9  and 
5.10,  were  10  to  46  dB  lower  than  with  the  orifice  exposed.   From 
these  figures,  the  major  contribution  from  the  orifice/louver  sources 
was  in  the  range  of  65  to  1000  Hz.   Background  noise  was  within  15  dB 
of  the  orifice/louver  noise  below  50  Hz  above  1000  Hz. 
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Figure  5.8:  Signal  to  noise  for  orifice  02  at  U  =  40  m/s. 
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Figure  5.10  Signal  to  noise  for  louver  11  at  U  =  50  m/s. 
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5.5   "Open  Window"  Contribution 

Previous  work  by  Splnka  [15]  determined  the  acoustic  pressure 
field  within  the  test  section  at  various  free  stream  velocities.   To 
determine  the  contribution  from  the  open  window  affect,  a  similar 
acoustic  pressure  field  was  generated  within  the  duct  using  a  speaker 
driven  a  B&K  random  noise  generator,  type  1402,  and  a  B&K  band  pass 
filter  set,  type  1612.   The  resulting  power  which  radiated  through  the 
orifice  was  measured  as  described  in  section  2.5  and  compared  to  the 
power  radiated  by  the  orifice  source.  Table  5.4  lists  the  resulting 
power  levels  for  both  cases  at  U  -  50m/s.   The  conclusion  was  that  the 
open  window  affect  results  in  levels  which  were  typically  23dB  lower 
than  the  levels  due  to  the  orifice  source. 


|  1/3  08  center 

Open  window 

Orifice  0-1 

frequency 

power  in  1/3  08  (2) 

power  in  1/3  08 

(3) 

(3-2) 

400 

53.8 

77.6 

23.8 

500 

48.2 

72.7 

24.5    | 

630 

45.9 

69.2 

23.3    | 

800 

43.7 

67.9 

24.2    | 

1000 

46.2 

66.4 

20.2    | 

1250 

41.0 

64.5 

23.5 

Table  5.4:  "Open  window"  power  levels  versus  orifice  0  -1  power  levels  at  U  =  50m/s 
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5.6     Orifice  Source  Power  Spectra  and  Total   Power 

The  measured  sound  power   spectra  versus    the   free   stream  velocity, 
U-20, 30,40   and   50  m/s ,    for  0-1,    0-2,    0-3,    L-l   and  L-2   are   shown  in 
Figures   5-11,    5-lla,    5-llb,    5-llc,    5-lld.      The   total   sound  power  radi 
ated  by  each  orifice/louver   source   at  the   above   speeds   is   given   in 
Table   5.5. 
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Figure  5.11     Orifice  01  sound  power  spectral   level  versus  free  stream  velocity 
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Figure  5.11(a):  Orifice  02  sound  power  spectral   levels  versus  freestream  velocity 
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Figure  5.11(b):  Orifice  03  sound  power  spectral    levels  versus  freestream  velocity 
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Figure  5.11(c):  Louver  L1  sound  power  spectra  levels  versus  freest  ream  velocity 
60. 

70.  - 

60.  " 
50.  ~ 


40.  ~ 


30. 


20.  ~ 


10. 


1000. 
FREQUENCY  (HZ) 


2000. 


Figure  5.11(d):  Louver  L2  sound  power  spectral  levels  versus  freestream  velocity 
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|    Orifice/      | 

Free  Stream 

Total  Sound  Power 

Louver 

Velocity  U(m/8)   | 

(dB  re  1  E- 12  watts)  | 

j     0-1 

20         | 

63.3 

30 

73.3         | 

40 

79.6         | 

SO 

84.0         | 

0-2 

20 

62.7 

30 

73.2         | 

40 

78.4 

50 

83.0         | 

0-3 

20 

62.9 

30 

72.2         | 

40 

76.3 

50 

79.7         | 

L-1 

20 

74.5 

30 

82.5         | 

40 

89.2 

50 

94.5 

L-2 

20 

64.9 

30 

73.6 

40 

80.1 

50 

85.3 

0-2A 

20 

62.0         | 

30 

71.8 

40 

77.2         | 

50 

82.2 

I-1A 

20 

58.1         | 

30 

68.4 

|         40 

73.3         | 

50 

77.7 

L-2A 

20 

55.3         | 

30 

65.0         | 

40 

69.7          | 

50 

76.3 

Table  5.5:     Orifice/ louver  total  sound  power  levels  versus  freestream  velocity,Ug 
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Several  observations  are  made  from  reviewing  the  above  data. 
5.6.1  The  effect  of  louverlng  on  sound  power 


(a)  louver,  L-l,  radiated  the  most  acoustic  power,  averaging  lOdB 
more  at  all  speeds  than  the  single  orifice  though  the  total 
open  areas  were  the  same  (see  Figures  5.12  and  5.12a). 

(b)  louver,  LI,  showed  strong  resonant  peaks  and  had  the  narrow- 
est half -power  bandwidth  of  all  the  samples  tested. 

(c)  louver,  LI,  radiated  significantly  more  power  than  louver  L2 
in  the  lower  frequency  bands  but  approximately  equal  power  in 
the  middle  to  high  frequency  bands  (see  Figures  5.13  and 
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2008. 


(d)  louver,  L2 ,  had  the  broadest  spectra  of  all  the  samples 
tested  with  equal  open  area. 

(e)  louver,  L2's  total  radiated  power  was  within  2  dB  of  orifice 
02 's  power  in  all  speed  ranges  but  the  spectral  shape  was 
significantly  broader.   As  shown  in  Figures  5.15  and  5.16, 
louver  L2  radiated  significantly  more  power  than  02  in  the 
middle  to  higher  frequency  bands  (10  dB  more  on  average)  but 
slightly  less  (approximately  1-2  dB)  power  in  the  lower 
frequency  bands . 
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Figure  5.15:  Effect  of   louven'ng  on  sound  power,   L2  versus  02  at  U  =  30m/s 
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Figure  5.16:   Effect  of   louven'ng  on  sound  power,   L2  versus  02  at  U  =  50  m/s 
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Thus  the  effect  of  louvering,  simply  by  dividing  the  total  open  area 

into  more  numerous,  smaller  areas  does  not  by  itself  reduce  the  total 

radiated  flow  noise  and  will  most  probably  result  in  increased  radi- 
ated sound  power  in  the  higher  frequency  bands. 


5.6.2  Orifice  throat  thickness  effect  on  radiated  power 

Figure  5.17  shows  the  orifice  throat  depth  effect  on  sound  power 
by  comparing  orifice  01  to  orifice  02.  Orifice  02  radiated  more  power 
(approximately  3dB)  in  the  higher  frequency  bands  but  maintained 
approximately  1-2  dB  less  power  in  the  more  dominant  lower  frequency 
bands . 
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Figure  5.17:  Orifice  depth  effect  on  sound  power,  01  versus  02  at  U  =  30  and  50  m/s 
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Thus  the  throat  thickness  was  more  predominant  in  effecting  the 
spectral  shape  than  reducing  the  total  radiated  power.   The  effect  of 
reducing  the  orifice  throat  thickness  appears  to  broaden  the  spectral 
shape  while  maintaining  the  same  total  radiated  sound  power. 

5.6.3  Streamwise  dimension,  b,  effect  on  sound  power 

The  spectral  shape  of  orifice  03  had  the  widest  half-power  band- 
width.  The  effect  of  reducing  the  streamwise  dimension,  b,  was  to 
broaden  considerably  the  power  spectra  while  reducing  only  moderately 
the  total  sound  power.   From  section  3,  the  controlling  shape  factor 
was  (l-cos(Kb)).  Thus  the  predicted  result  from  reducing  b  by  four 
fold  would  be  a  subsequent  increase  by  four  fold  in  the  power  spectra 
bandwidth.   Reviewing  again  the  power  spectra  for  orifices  02  and  03 
(Figures  5.11a  and  5.11b),  we  see  that  indeed  the  power  bandwidth  for 
03  is  approximately  four  times  that  of  02  (defining  the  bandwidth  to 
be  between  the  one-quarter  peak  spectral  power  points). 

5.7  Speed/Power  Laws 

The  speed/power  law  is  of  the  form: 

P  -  C  Un 
where  C   is  a  value  which  is  independent  of  speed.   The  data  in  Table 
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5.5  was  utilized  to  determine  the  total  radiated  power  dependence  on 
speed.  The  equations  used  to  determine  the  value  of  the  speed  expo- 
nent, n,  were: 

P3/P2  -  (u3/u2)n 

10  log  (P3/Pref)  -  10  log(P2/Pref)  "  n  10  log(U3/U2) 

Lp3-  Lp2"  *  log(U3/U2) 

n  -  (Lp3-  Lp^/(10  log  (U3/U2)) 

where  the  subscripts  3  and  2  correspond  to  30  and  20  m/s.   Similar 
equations  hold  for  U4/U2,  U5/U2,  U4/U3,  U5/U3  and  U5/U4.   Table  5.6 
lists  the  calculated  value  for  n  for  each  sample  tested.   The  values 
for  n  varied  from  3.3  to  6.8  with  a  mean  value  of  4.9  and  a  standard 
deviation  of  0.3.   Excluding  the  extremely  low  values  of  n  recorded 
for  03  at  the  higher  free  stream  velocities,  the  mean  value  of  n  was 
5.0  with  a  standard  deviation  of  0.2. 

The  speed/power  relationship  for  orifices  01,  02  and  louvers  LI 
and  L2  is  plotted  in  Figure  5.18.   Also  shown  is  the  least  squares  fit 
to  the  data  points  satisfying  the  equation: 

(LPi  -  LPj)  -  n  10  logOJi/Uj)  (5.2) 

From  the  fit,  the  best  value  for  n  is  4.97. 
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Or  if  ice/Louver 

U3/U2 

Vu2 

Vu2 

U4/Uj 

U5/U3 

u5/u4 

|           0-1 

5.7 

5.4 

5.0 

5.2 

4.8 

4.5                  | 

|           0-2 

6.0 

5.2 

4.2 

5.1 

4.4 

4.7                  | 

|           0-3 

5.3 

4.5 

3.3 

4.2 

3.4 

3.5                  | 

I           L"1 

4.6 

4.9 

5.4 

5.0 

5.4 

5.5                 | 

I           L'2 

4.9 

5.2 

5.5 

5.1 

5.3 

5.0                 | 

0-2A 

5.6 

5.1 

4.3 

5.1 

4.7 

5.2                 | 

L-M 

5.9 

5.1 

3.9 

4.9 

4.2 

4.5                | 

|           L-2A 

5.5 

4.8 

3.8 

5.3 

5.1 

6.8                 | 

average 

5.4 

5.0 

4.5 

5.0 

4.7 

5.0 

std.  dev. 

0.5 

0.3 

0.6 

0.3 

0.6 

0.9 

overall  average 

4.9 

overall 

standard  deviation 

0.3 

Table  5.6   Speed/Power  relationship  (P  s  CU°),  values  for  the  speed  exponent,  n 
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Thus  the  concluding  dependence  is: 

Total  radiated  power  increases  with  the  5th  power  of  the  free 
stream  velocity 

5.8  Directivity  of  orifice  01  and  louver  L2A 

Directivity  measurements  from  orifice  01  and  louver  L2A  were 
taken  at  U  -  25m/s.   Figure  5.19  details  the  coordinate  system  used 
for  the  measurements .   The  measurements  indicated  the  source  was 
nearly  omni -directional.   Power  levels  in  the  XZ  (azimuth)  plane  dev- 
iated by  0.9  dB  from  the  maximum  level.  In  the  YZ  (polar)  plane,  power 
levels  deviated  by  2.6  dB  from  the  maximum  level.  Off  the  major  axis, 
power  levels  deviate  by  3.6  dB  from  the  maximum  level  directly  above 
the  orifice  source.   Table  5.7  lists  the  measured  sound  intensity  in 
dB  re  IE- 12  watt/m^  at  a  radial  distance  of  r  -  27  inches  and  the 
various  polar  and  azimuth  angles.   Figure  5.20  is  the  azimuth  direc- 
tivity plot  for  orifice  01  and  louver  L2A  and  clearly  supports  the 
conclusion  of  a  nearly  omni-directional  source. 
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Figure  5.19  Coordinate  system  for  directivity  measurements.  Reference  is  the 
aperture  center.  Polar  angle,  9,   measured  from  Z  axis,  azimuthal 
angle,  4,   measured  from  -x  axis.  Measurements  described  by  spheri- 
cal coordinates,  r,9,4. 


Total 

Intensity 

Polar  angle 

Azimuthal 
Orifice 

angle 
0-1 

dB  re 

lE-^watt/m2  j 

0 

90 

57.3 

0 

45 

56.9 

0 

15 

56.4 

0 

135 

57.1 

0 

165 

56.8 

0 

105 

57.2 

0 

75 

57.0 

I     -85 

90 

54.9 

•45 

90 

56.5 

0 

90 

57.5 

45 

90 

57.0 

85 

90 

55.9 

-45 

30 

54.7 

•45 

25 

53.9 

•45 

35 

55.0 

Table  5.7  Directivity  measurement  for  orifice  0-1  and  louver  L2A 
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Total  Intensity 

Polar  angle 

Azimuthal  angle 
Louver  L2A 

d8  re  1E-12watt/m2  | 

0 

160 

66.2       j 

0 

145 

66.1 

0 

130 

66.0 

0 

115 

66.0 

0 

90 

65.8       | 

0 

75 

65.7       | 

0 

60 

65.7 

0 

45 

65.8       | 

0 

30 

65.2       | 

0 

15 

65.2       | 

•90 

90 

63.4       | 

•75 

90 

63.9 

-60 

90 

65.0 

•45 

90 

65.4 

•30 

90 

65.9 

-15 

90 

66.1 

0 

90 

66.2 

15 

90 

66.4 

30 

90 

66.5 

45 

90 

66.6 

60 

90 

66.4 

75 

90 

66.3 

90 

90 

65.6 

•45 

45 

65.3 

|     -45 

30 

64.2 

-45 

90 

65.3       | 

-45 

135 

66.0       | 

Table  5.7  Directivity  measurement  for  orifice  0-1  and  louver  L2A  (Continued) 
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Figure  5.20  Orifice  01  directivity  in  the  XZ  (azimuth)  plane  and  louver  L2A 
directivity  in  the  YZ  (polar)  plane 
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5.9   Edge  Angle  Effect  on  Radiated  Sound  Power 

5.9.1  Introduction 

Earlier  It  was  shown  that  Indeed  the  orifice  source  provided  both 
the  spectral  character  and  the  predominant  amount  of  the  total  power. 
The  strength  of  this  source  is  due  in  part  to  the  eddy  interaction 
with  the  downstream  edge  of  the  orifice.   To  determine  the  relative 
importance  of  this  edge  interaction  compared  to  the  other  significant 
parameters  identified  in  Section  3,  three  additional  orifice/louver 
samples  were  tested.   They  were  orifice  0-2A  and  louvers  L-1A  and 
L-2A.   The  details  of  these  samples  are  described  In  Section  4.   The 
only  difference  between  these  samples  and  the  previously  tested 
samples  was  the  leading  and  trailing  edges  were  slanted  away  from  the 
flow  by  40  degrees.   This  provided  both  a  blunt  trailing  edge  where 
the  eddies  strike  and  a  sharpened  leading  edge  (essentially  reducing 
the  upstream  aperture  to  1/32",  the  radius  of  the  edge). 

5.8.2  Angled  edge  orifice/louver  power  spectra 

The  measured  sound  power  spectra  for  the  above  samples  are  shown 
in  Figures  5.21,  5.22  and  5.23.   The  total  power  radiated  by  each 
sample  at  the  various  freestream  velocities  is  given  in  Table  5.5. 
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Figure  5.21:     Orifice  02A  power  spectra  level  versus  freestreaa  velocity 
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Figure  5.22:     Louver  L1A  power  spectral   levels  versus  freestream  velocity. 
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Figure  5.23:  Louver  L2A  power  spectral  levels  versus  freestrean  velocity. 

Several  observations  are  worth  noting: 

(a)  the  effect  of  angling  the  edges  is  more  predominant  for  the 
louvers  than  the  single  orifice. 


(b)  the  effect  on  louver  L-1A  is  a  14-16  dB  reduction  in  total 

radiated  sound  power.   The  large  resonant  peaks  which  occurred 
in  the  spectra  for  L-l  are  not  present  in  the  spectra  for 
L-1A.   Figures  5.24  and  5.25  compare  the  two  louvers  at  U=20 
and  50  m/s.   Not  only  are  the  resonant  peaks  eliminated  by 
angling  the  edges  but  also  the  radiated  power  is  reduced  by 
5-7  dB  in  the  higher  frequency  bands. 
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Figure  5.24:  Edge  angle  effect  on  sound  power  spectral  level,  LI  versus  LU 
at  U  *  20  m/s 
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Figure  5.25:   Edge  angle  effect  on  power  spectral    level,    L1  versus  L1A  at 
U  =  50  m/s 
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(c)  Figure  5.26  compares  the  power  spectral  levels  of  L-2  and 
L-2A.   The  result  for  this  louver  is  a  8.6  to  10.4  dB  reduc- 
tion in  sound  power  equally  divided  across  the  frequency 
bands . 

(d)  for  the  single  orifice,  0-2A,  the  effect  of  angling  the  edges 
resulted  in  only  a  0.7  to  1.4  dB  reduction  in  sound  power. 
Figures  5.27  and  5.28  compare  the  power  spectra  for  orifice 
0-2  and  0-2A  at  U-20  and  50  m/s .   Here  again  the  effect  of 
angling  the  edges  resulted  in  a  reduction  in  sound  power  in 
the  higher  frequency  bands  but  in  this  case,  essentially  no 
change  in  the  lower  frequency  bands . 
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Figure  5.26:  Edge  angle  effect  on  sound  power  spectral  level,  12  versus  L2A  at 
U  =  50  m/s. 
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Figure  5.27:     Edge  angle  effect  on  sound  power  spectral   level, 
02  versus  02A  at  U  =  50  m/i. 
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Figure  5.28:  Edge  angle  effect  on  sound  power  spectral  level, 
02  versus  02A  at  U  =  20  m/s. 
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5.9.3   Combined  effect  of  louverlng  and  angling  edges 

Taken  together,  the  effect  of  louvering  and  angling  edges  results 
in  a  reduction  of  3.4  to  7.5  dB  in  radiated  sound  compared  to  the 
single  angled  orifice  (both  of  equal  total  open  area) .   This  reduction 
most  probably  is  caused  by  two  factors: 

(a)  the  louvers  shed  eddies  which  after  convecting  downstream  to 
the  trailing  edge  are  the  same  order  in  size  as  the  edge 
thickness.   This  results  in  a  more  beneficial  interaction  when 
the  eddies  impact  the  angled  edge.   The  orifice,  however,  with 
a  four  times  larger  aperture,  sheds  eddies  which  are  much 
larger  in  size  upon  impact  with  the  trailing  edge  and  are 
therefore  less  affected  by  the  thin  angled  edge. 

(b)  louvering,  by  reducing  the  transverse  dimension,  b,  interferes 
with  the  longitudinal  correlation  of  the  shed  eddies.   This 
interference  will  be  more  predominant  in  the  lower  frequencies 
where  the  longitudinal  correlation  length  is  larger  than  the 
aperture  transverse  dimension,  b.   As  shown  in  Figures  5.29, 
5.30,  5.31  and  5.32,  the  reduction  in  sound  power  is  predom- 
inantly in  the  lower  frequency  bands.   There  is,  however,  more 
radiated  power  in  the  angled  louver  in  the  higher  bands  by 
approximately  1-3  dB. 
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Figure  5.29:     Effect  of  Louvering  on  sound  power  spectra  level,   L2A 
versus  02A  at  U  *  30  m/s 
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Figure  5.30:  Effect  of  louvering  on  sound  power  spectral  level,  L1A 
versus  02A  at  U  =  30  m/s 
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Figure  5.31:     Effect  of   louvering  on  sound  power  spectral    level,   12A  versus  02A 
at  U  =  SO  m/s 
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Figure  5.32:     Effect  of    louvering  on  sound  power  spectral    level,   L1A 
versus  02A  at  U  =  50  m/s. 
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Thus  one  can  conclude  that  effective  louvering  does  reduce  the 
total  radiated  power  by  reducing  single  mode  resonance  and  by  inter- 
fering with  the  lower  frequency,  longitudinal  eddy  correlation.   This 
results  in  effectively  broadening  the  power  spectrum. 

5. 9. A  Louver  geometry  effect  on  sound  power 

As  with  the  straight  walled  louvers,  the  multiple  circular  aper- 
ture louver  was  quieter  than  the  multiple  rectangular  aperture  louver 
by  1.4  to  3.6  dB  in  total  radiated  power.   As  shown  in  Figures  5.33 
and  5.34,  louver  L-2A  had  approximately  the  same  spectral  shape  as 
louver  L-lA  but  with  considerably  less  power  in  the  low  frequency 
bands.   Again  this  reduction  can  be  attributed  to  the  additional 
interference  of  the  lateral  correlation  due  to  the  reduced  dimension 
in  the  aperture's  height.   To  test  this  conclusion,  one  can  calculate 
the  frequency  where  the  lateral  correlation  length,  i,  equals  the 
aperture  diameter,  d  for  louver  L-2A  and  at  freestream  velocity  of  50 
m/s.   From  Figure  3.4,  at  B(K(r3>)-0.4,  the  value  of  Kr3  equals  1.2. 
Setting  r3  equal  to  0.01  m  (27/64")  and  U-50  m/s  gives  a  limiting 
frequency  of  450  Hz.   Thus  one  would  predict  diminishing  reductions  in 
sound  power  for  frequencies  above  450  Hz.   Reference  to  Figure  5.34 
shows  that,  indeed,  at  frequencies  above  450  Hz  the  effect  of  lateral 
correlation  interference  on  reducing  sound  power  is  diminishing. 
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Figure  5.33:  Louver  geometry  effect  on  sound  power  spectral  level, 
L1A  versus  L2A  at  U  *  30  m/s. 
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Figure  5.34:     Louver  geometry  effect  on  sound  power  spectral    level, 
L1A  versus  L2A  at  U  =  50  m/s. 
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Figure  5.35:  Turbulent  intensity  spectral  levels  for  louver  LI  versus  freestream 
velocity. 
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5.10  Describing  the  Velocity  Field  Surrounding  the  Samples 


In  Section  3,  the  radiated  sound  power  was  shown  to  be  directly 
related  to  the  normal  component  of  the  turbulent  intensity.   Several 
measurements  of  the  turbulent  intensity  spectrum  were  taken  as 
described  in  Section  4.4.   The  assumption  was  made  that  the  turbulent 
intensity  in  the  streamwise  and  normal  component  directions  were 
equal.   Thus  3  dB  was  subtracted  from  the  measured  intensity  spectral 
levels  to  give  only  the  normal  component  of  the  turbulent  fluctua- 
tions. The  following  observations  are  made: 
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(a)  the  turbulent  Intensity  measured  2  tenths  of  an  Inch  Into  the 
boundary  layer  was  relatively  unaffected  by  the  probe  loca- 
tion in  the  transverse  direction  and  the  various  orifice 
geometries.   Figure  5.35  is  the  measured  normal  component  of 
the  turbulent  intensity  for  louver  L-l,  with  the  probe  at  the 
center  of  the  upstream  aperture.   These  spectral  levels  were 
used  in  the  following  section  to  collapse  the  sound  power 
spectra. 

(b)  there  were  significant  reductions  in  turbulent  intensity 

(approximately  3  dB  )  between  louvers  L-l  and  L-lA  when  meas- 
ured at  the  leading  edge  of  the  leading  aperture.  This  reduc- 
tion in  intensity  supports  the  measured  reduction  in  total 
sound  power  and  also  supports  the  conclusions  that  the  sound 
power  is  related  to  the  driving  aerodynamic  field. 

Despite  a  fairly  rigorous  attempt  to  maintain  the  block-house  air 
tight,  a  finite  and  measurable  pressure  differential  existed  between 
the  blockhouse  and  the  duct  static  pressure.  This  pressure  differen- 
tial resulted  in  a  small  amplitude,  mean  flow  entering  the  ori- 
fice/louver from  the  closed  acoustic  space.   The  magnitude  of  this 
transverse  mean  flow  was  greatest  in  the  louver  samples  and  measured 
on  average,  fourteen  percent  of  the  freestream  velocity.   The  exact 
contribution  of  this  mean  flow  to  the  radiated  sound  power  in  not 
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known  but  It  is  estimated  to  be  of  minimal  impact  since  the  turbulent 
intensity  of  the  entering  flow  measured  several  orders  of  magnitude 
less  than  the  normal  component  of  the  freestream  turbulence. 

5.11  Verifying  the  Scaling  Laws 

In  Section  3,  an  equivalent  source  model  was  developed  for  rectan- 
gular shaped  aperture  orifices  and  louvers.   To  ascertain  the  validity 
of  this  model,  the  sound  power  spectra  for  orifices  02,  03  and  the 
louver  LI  were  collapsed  by  normalizing  the  spectral  levels  by  the 
magnitude  of  the  equivalent  source  model.   The  resulting  collapsed 
spectra  are  shown  in  Figures  5-36,  5-37  and  5-38.   It  is  evident  by 
the  collapse  that  the  equivalent  source  model  adequately  describes  the 
controlling  parameters  for  orifice  and  louver  flow  noise.   The  predic- 
tion of  the  actual  magnitude  of  the  sound  power  is  fairly  accurate  for 
the  samples  with  the  smaller  apertures  but  is  too  high  by  about  a  fac- 
tor of  100  for  the  larger  aperture  openings  (01  and  02).   The  probable 
reasons  for  this  are  (1)  overestimating  the  effective  aperture  width, 
b,  (2)  the  estimation  for  the  lateral  correlation  lengths  and  (2)  the 
assumption  of  equal  components  of  turbulent  intensity.   As  mentioned 
in  Section  3,  for  orifice  0-2,  the  longitudinal  correlation  amplitude 
function,  A(K(c^)  does  attenuate  to  a  value  of  approximately  0.6  for 
the  slower  speeds  indicating  that  some  eddy  decay  occurred  in  the 
aperture.   For  this  case,  one  could  assume  a  smaller  transverse 
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length,  b' ,  over  which  the  eddy  decay  was  minimal  and  use  this  length 
scale  for  evaluating  the  source  strength.   Defining  b'  at  the  point 
where  the  magnitude  of  the  amplitude  function  equaled  0.95  results  in 
an  effective  aperture  length  of  0.5  inches  versus  the  actual  dimension 
of  b-1.67  inches  (for  orifice  0-2  at  U-20m/s  and  f-250  Hz).   This 
results  in  a  reduction  of  4.2  dB  in  the  predicted  spectral  level  and 
could  account  for  some  of  the  error  in  predicting  the  larger  aperture 
radiated  sound  power. 
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Figure  5.36:     Scaling  law  verification  for  orifice  02 
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Figure  5.37:     Scaling  taw  verification  for  orifice  03 
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Figure  5.38:     Scaling  law  verification  for  louver  L1 
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6.   CONCLUSIONS 

The  acoustic  source  created  by  turbulent  flow  over  orifices  and 
louvers  was  investigated  by  measuring  the  radiated  sound  Intensity, 
utilizing  an  intensity  probe,  at  various  wind  speeds  In  a  semi- 
anechoic  wind  tunnel.   From  the  intensity  spectra,  the  radiated  power 
and  the  directivity  of  the  source  was  determined.   The  goal  of  this 
research  was  to  describe  the  source  and  identify  the  controlling  par- 
ameters which  affect  its  spectral  shape,  radiated  power  and  directi- 
vity.  The  orifice  and  louvers  investigated  have  geometries  similar  to 
those  used  in  industry  in  both  air  and  water  environments.   The 
effects  of  various  leading/trailing  edge  geometries  were  also  investi- 
gated. 

The  source  of  noise  has  been  modelled  on  the  basis  of  the  theoret- 
ical analysis  by  Ffowcs  Williams,  Nelson  and  Corcos.   The  equivalent 
source  of  noise  was  modelled  as  two  equal  but  opposite  monopole 
sources  located  In  the  upper  and  lower  aperture  of  each  orifice.   The 
strength  of  each  source  is  determined  by  the  mass  flow  driven  through 
the  apertures  by  the  normal  component  of  the  turbulent  fluctuations. 
The  resulting  radiated  sound  field,  both  into  a  closed  acoustic  space 
and  radiating  away  from  the  orifice,  is  dependent  on  the  acoustical 
properties  of  the  louver. 
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To  ascertain  the  validity  of  the  source  model,  two  orifice  and  one 
louver  power  spectra  were  collapsed  by  normalizing  the  spectral  levels 
by  the  magnitude  of  the  equivalent  source  model.   It  is  evident  by  the 
collapse  that  the  equivalent  source  model  adequately  describes  the 
controlling  parameters  for  orifice  and  louver  flow  noise.   The  con- 
trolling parameters  for  the  geometries  tested  were: 

(1)  Strouhal  number  based  on  the  lateral  correlation  length,  i, 
and  aperture  transverse  dimension,  b 

(2)  the  normal  component  of  the  Turbulent  Intensity, 

(3)  the  fifth  power  of  the  freestream  velocity  and 

(4)  the  shape  factor  (l-cos(Kb))  where  K  is  the  convective 
wavenumber,  K  -  w/Uc,  and  Uc  is  the  eddy  convection  velocity. 

As  a  result  of  this  research,  several  important  conclusions  are 
made : 

(a)  sound  intensity  measurements,  using  an  intensity  probe,  can 
successfully  describe  an  acoustic  source,  even  though 
measurements  are  taken  within  the  source's  near  field. 

(b)  the  errors  associated  with  measuring  intensity,  using  an 
intensity  probe,  can  be  as  small  as  the  inherent  statistical 
deviation  of  the  random  fluctuations. 
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(c)  the  equivalent  source  model  developed  to  describe  the  behav- 
ior of  orifice  and  louver  flow  noise  successfully  describes 
the  controlling  parameters  and  provides  valuable  insight  for 
predicting  the  actual  magnitudes  of  the  radiated  sound  spec- 
trum. 

(d)  the  1/3  octave  band  containing  the  most  acoustic  power  from 
orifice  and  louver  flow  noise  can  be  predicted  from  the  non- 
dimensional  Strouhal  number  based  on  the  aperture's  trans- 
verse dimension  and  the  non-dimensional  length  scale,  S/b 
where  6   is  the  turbulent  boundary  layer  thickness. 

(e)  for  aerodynamic  flows,  the  orifice/louver  source  provides 
both  the  dominant  amount  of  total  sound  power  and  the  spec- 
tral character  of  the  radiated  noise.   Contributions  due  to 
structural  radiation  from  the  ducting  and  ducting  "open  win- 
dow" effects  are  not  significant. 

(f)  the  effect  of  louver ing,  simply  by  dividing  the  total  open 
area  into  more  numerous,  smaller  areas  does  not  by  itself 
reduce  the  total  radiated  flow  noise  and  will  most  probably 
result  in  increased  radiated  sound  power  in  the  higher 
frequency  bands . 

(g)  the  effect  of  reducing  the  aperture  throat  thickness  is  to 
broaden  the  spectral  shape  of  the  radiated  power  spectrum 
while  maintaining  approximately  the  same  total  radiated  sound 
power. 
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(h)   the  effect  of  reducing  the  aperture's  transverse  dimension 
can  be  adequately  described  by  the  controlling  shape  factor 
in  the  equivalent  source  model,  (l-cos(Kb)). 

(i)  directivity  measurements  show  a  nearly  omni-directional 
source  in  the  on-axis  planes  with  off-axis  power  levels 
reduced  by  2-3  dB. 

(j)   taken  together,  the  effect  of  louvering  and  angling  the  lead- 
ing and  trailing  edges  resulted  in  reductions  to  sound  power 
of  3.4  -7.5  dB  compared  to  the  equivalent  area,  single  aper- 
ture, angled  orifice. 

(k)   the  measured  reductions  in  sound  power  due  to  effective  lou- 
vering can  be  explained  by  (1)  reductions  in  the  normal  com- 
ponent of  the  turbulent  intensity  due  to  the  eddies  impacting 
a  blunt  trailing  edge  and  (2)  destructive  interference  with 
the  eddy  longitudinal  and  lateral  correlations. 

(1)   the  louvers  with  numerous  circular  apertures  and  angled  edges 
radiated  the  least  acoustic  power  with  a  very  broad,  flat 
power  spectrum. 
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APPENDIX  1:  MATCHED  ASYMPTOTIC  EXPANSION  SOLUTION  TO  A  RADIATING 
ORIFICE  IN  AN  INFINITE  RIGID  WALL 

The  two  dimensional  problem  to  be  solved  is  the  far  field  radiated 
pressure  at  x^,  created  by  a  driving  aerodynamic  field  at  y^,  within  a 
finite  size  aperture  in  an  infinite,  rigid  wall.   Using  the  principle 
of  reciprocity,  the  above  field  can  be  determined  by  first  solving  for 
the  sound  field  at  y^  created  by  a  distant  monopole  at  x^.   Then,  the 
desired  solution  is  obtained  by  simply  replacing  the  position  vectors 
of  the  far  field  monopole  and  the  aperture  source.   The  solution  to 
this  problem  will  be  analyzed  using  the  principles  of  Matched  Asymp- 
totic Expansions  (MAE) . 


The  coordinate  system  used  in  the  analysis  is  shown  below, 

% 


-b 


\ 


with  boundary  conditions,  3p/3y2  -  0  for  |y^|  >  b  and  y2  -  0. 


The  governing  equations  for  a  plane  wave  of  single  frequency,  u>, 
at  normal  incidence  are: 
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V2p  +  k2p  -  0  (Al.l) 

l(ky  -wt) 
p  -  Ae  (A1.2) 

where  k  is  the  acoustic  wavenuniber,  k  -  <*>/cq.   Also  define  the 
non-dimensional  terms,  r  -  wt,  y^  -  ky^,  y2  -  ky2 ,  p  -  p/A  and  e   -  kb 
for  the  parameters  in  the  outer  problem. 

Substituting  the  above  outer  parameters  into  equations  Al.l  and 
Al.2  gives  the  governing  equations  and  boundary  condition  for  the 
outer  problem 


V2p  +  p  -  0  (A1.3) 


Ky2-0 

P  -  e  .  (A1.4) 


To  lowest  order,  c-»0,  and  the  orifice  closes.   Thus  the  outer 
problem  becomes 

p  -  $lPl  +  $2P2  +  ***  as  €~*®  (Al.5) 

with  6i   -  1.   From  the  boundary  condition,  we  can  determine  p^  to  be 


iy2    -iy2 
Pi  ~  Pi  +  Pr  "  e    +  e    •  v2  <  °  (Al.6) 
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suppressing  the  e**r  factor,  where   p^,  pr  are  the  incidence  and 
reflected  waves,  respectively.   Note  that  the  boundary  condition  is 
satisfied 

dpi  iy2   -iy2 

-  i  (e   -e    )  -  0  for  y2  -  0  (Al.7) 

3y2 

Choose  also  the  condition  that  p^  ■  0  for  y2  >  0.   To  open  the 

orifice,  choose  the  inner  variables  y^  -  y\/t   and  y2  -  y2/« . 

In  inner  variables, 

p  -  exp(i(cy2-r)]  +  exp[ -i(«y2-r) ] 

-  exp(-ir){exp(ey2)  +  exp(-«y2)} 

2~2  2~2 

—         —  *  t\  —  J  n 

-  exp(-ir) {l+cy2  -  *  -  •••+!  -  cy2  -  *  +•••} 

2  2 

2"2 

_        e  y  _ 

-  2e"lf  {1-  2  +  G(«4))  as  <  -  0,y2  fixed  (y2  -  0)     (Al.8) 

2 

The  governing  equation  with  the  Inner  parameters  substituted  into 
equations  Al.3  becomes 

V2p  +  «2p  -  0  (A1.9) 

with  boundary  conditions  3p/3y2  -  0  for  y2  -  0  and  \y\\    >  1.   Let 
p  be  defined  by 

p  -  A0(Opi  +  A].(OP2  +   '"  (ALIO) 
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for  the  inner  expansion.   The  resulting  inner  problem  coordinate  sys 
tem  is  i  /2. 


-l 


Yi 


Since  V^p^  -  0,  solve  the  inner  problem  by  conformal  mapping,  p^  - 

A  Rew.   Using  a  solution  of  the  form,  iz  —  cosh  w, (suppressing  the 

e'ir  factor)  we  can  determine 

1     1-3 

w  -  in  2  z  +  iw/2  + +•••  as  |z|-«*>  (Al.ll) 

_2       _4 
2x2z   2x4x4z 


where  z  -  re 


LB 


Since  z  -  z/e  , 


w  -  in  2z  -  inc  +  iw/2  +  $(«2)  for  e   -  0,  z  fixed  (z-"»)   (Al.12) 


Therefore,  utilizing  equation  A. 11,  we  can  determine  Aq(€)pi  to  be 


A0(OPi  -  AA0{in  2r  -  ine  +  iw/2  +  Q(e2)) 


(A1.13) 
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Now  matching  the  leading  terms  in  the  outer  and  inner  solutions, 
equations  Al.8  and  Al.13,  leads  to  Aq  -  +  1/in  «  and  A  -  -  2.   Then 
$2(<)  ~  1/^n  e    in  the  outer  expansion.   P2  satisfies  the  Helmholtz 
equation  with  k  -  1  and  must  behave  like  In  r  for  r  -»  0,  c   -»  0  and  r 
fixed.   Clearly 

(1) 
p2  -  C  H0  (r)  (A1.14) 

satisfies  this  requirement  and  also  satisfies  the  rigid  wall  condi- 
tion. As  r  -*  0, 

(1)      21  2i 

Hq  (r)  -  —  in  r  +  1  +  —  (7  -  in2)+««» 

n  n 

where  7  is  defined  as  Euler's  constant.   Thus  to  match  with  the  inner 
solution,  A  -  -  2  -  C  •  2i/»,  C  must  be  equal  to  i*.   Therefore  the 
final  outer  expansion  becomes 


,    v2    "v2   i*He(1)(r) 
p  -  e"lr  {e   +  e    + +  •••  }  (Al.15) 

ine 


where  the  first  term  in  brackets  is  the  incident  wave,  the  second  term 
is  the  reflected  wave  and  the  lowest  order  orifice  effect  is  the  last 
term  which  describes  a  monopole  type  radiation.   Substituting  for  the 
non-dimensional  terms  and  applying  the  principle  of  reciprocity,  one 
obtains  the  radiated  sound  from  a  2  dimensional,  single  slot,  line 
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source 


kx2    -kx2   iwHJ^Cr)   -Iwt 

P  -  A{e    +  e     +  )e  (Al.16) 

inkb 


(1)  -    I r*  i(r-w/4) 

As  r  -**>,  H   (r)  -♦  j2/*r  e         so  the  for  field  radiated 


sound  for  the  2D,  single  slot,  line  source  becomes 


kx2     -kx2      i =n  i(r-w/4)   -iwt 

P  -  A{e     +  e      +  iwj2/wr  e         )e  (A1.17) 

r-«o 


Equation  Al.17  provides  the  result  that  the  far  field  radiated 
pressure  for  a  2D,  single  slot  source  to  lowest  order  is  composed  of: 

(1)  the  incident  and  rigid  wall  reflection 

(2)  2D  line  source  having  a  1/r  dependence  (cylindrical  spread- 
ing). 

Thus  one  can  infer  that  it  is  plausible  for  the  3  dimensional, 
single  orifice  source  to  have  a  reflection  coefficient,  |R|  ~   1  for 
values  of  b/A  «  1. 
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APPENDIX  2:   MEASURING  TURBULENT  INTENSITY 

In  this  appendix,  the  transfer  function  to  relate  the  mean  and 
fluctuating  component  of  the  hot  wire  anemometer  output  voltage  to  the 
mean  and  fluctuating  component  of  the  particle  velocity  is  developed. 

The  basic  relationship  between  output  voltage,  E,  and  particle 
velocity,  U,  is  defined  by  King's  Law 

E  -  A  +  BUn  (A2.1) 

Now  define  the  instantaneous  voltage  and  particle  velocity  by 

E  -  E0  +  e',  |e'|  «  |E0|  (A2.2) 

U  -  U0  +  u\  |u'|  «  |U0|  (A2.3) 

where  Eq  is  the  mean  of  E,  Uq  is  the  mean  of  U,  e'  and  u'  are  the 
fluctuating  components  whose  expectations  are  zero.   Substituting 
equations  A2 . 2  and  A2  .  3  into  equation  A2 . 1 ,  expanding  and  dropping  the 
higher  order  terms  of  u'  and  e'  gives 


E2  +  2E0e'  -  A  +  BUn  (1  +  [nu'/U0])  (A2.4) 


Taking  the  expectation  of  equation  A2.4  gives 
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E2  -  A  +  BUn  (A2.5) 

0  0 


Subtracting  equation  A2 . 5  from  equation  A2.4  gives  the  desired 
relationship  between  the  fluctuating  components  of  the  anemometer 
output  voltage,  e',  and  the  particle  velocity,  u' , 

1-n 
u»  -  ([2E0U0   ]/Bn)e'  (A2.6) 

Thus,  defining  the  term  in  parentheses  as  H,  the  resulting  rela- 
tionship between  the  spectral  densities  of  the  fluctuating  velocity, 
Su(f),  and  the  output  voltage,  Sj?(f),  is 

Su(f)  -  |H|2  SE(f)  (A2.7) 
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APPENDIX  3:   MATHEMATICAL  STEPS  IN  DETERMINING  SQ(yk,u) 

From  equation  3.7,  the  cross  spectral  density  is  assumed  to  be  of 
the  form 

Sw((k,u)  -  Sv(w)A(K€1)B(K€3)  •  exp  (iK^)  (A3.1) 

As  described  in  Section  3.2,  the  normalized  longitudinal  cross - 
spectra  amplitude  function,  A(K«i),  is  identically  equal  to  one  for 
all  samples  except  01  and  02  at  the  lower  freestream  velocities,  where 
at  a  minimum  it  equals  0.6.   So  assuming  no  decay  in  the  eddy  struc- 
ture in  the  y^  direction  (i.e.,  A(K«^)  equals  one)  and  Sv(yk,w)  to  be 
uniform  over  the  aperture  surface,  the  source  strength  spectral  den- 
sity becomes 

SQ(yK.">.  "  A>Vsv(yk,W)/jB(Ke3)exp(iK€l)tf(€k)dA(€k)       (A3. 2) 

00 

where  c^,  €3  are  the  position  coordinates  as  shown  below. 
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Substituting  for  0(«k)  we  have 


ike, 

SQ<yk,w)  -  P1"2   Sv(yk,u>)  /  /  B(K<3)  e    (L-|«3|)(b-|e1|)   (A3. 2a) 


•  dA  (ek). 

Since  the  lateral  correlation  1,  is  typically  one- tenth  or  less 
the  length  of  the  aperture,  the  integral  in  the  €3  direction  can  be 
evaluated,  independent  of  the  actual  limits  of  the  integrand,  as 


I 


B(K€3)(L-|«3|)d€3  -  Li  -  0(i2) 


(A3. 3) 


where 


I 


B(K«3)dc3  -  i  (A3. 3a) 


Evaluating  the  integral  over  the  c^  directions  gives: 


I 


iKc.  2 


e   1  (b-\ei\)6el   -  —  (l-cos(Kb))  (A3. 4) 


K2 


Therefore  the  final  result,  neglecting  the  i2  term  is: 

2 
Sn(yk,W)  "  pVsv(yk,a))Li(_)  (l-cos(Kb))  (A3. 5) 

K2 
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The  remaining  terra  needed  to  fully  describe  the  magnitude  of  the 
source  strength  spectral  density  is  i.  The  lateral  correlation  length, 
i,  is  determined  at  the  point  where  the  amplitude  function, 
B(K«3>  equals  0.4  and  the  frequency  equals  the  center  frequency,  f c , 
as  defined  in  Section  5.3. 
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